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The heavy fourth generation of quarks that have sufficiently small mixing with the 
three known SM families form hadrons. In the present work, we calculate the masses 
and decay constants of mesons containing either both quarks from the fourth gener- 
ation or one from fourth family and the other from known third family SM quarks in 
the framework of the QCD sum rules. In the calculations, we take into account two 
gluon condensate diagrams as nonperturbative contributions. The obtained results 
reduce to the known masses and decay constants of the bb and cc quarkonia when 
the fourth family quark is replaced by the bottom or charm quark. 
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I. INTRODUCTION 

In the standard model (SM), we have three generation of quarks experimentally observed. 
Among these quarks, the top (t) quark does not form bound states (hadrons) as a conse- 
quence of the high value of its mass. The top quark immediately decays to the bottom 
quark giving a W boson and this transition has full strength. The number of quark and 
lepton generations is one of the mysteries of nature and can not be addressed by the SM. 
There are flavor democracy arguments that predict the existence of the fourth generation 
of quarks [IHiJ]. It is expected that the masses of the fourth generation quarks are in the 
interval (300 — 700) GeV in which the upper limit coincides with the one obtained from 
partial- wave unitarity at high energies jBf. Within the flavor democracy approach, the Dirac 
masses of the fourth family fermions are almost equal, whereas masses of the first three 
family fermions as well as the CKM and PMNS mixings are obtained via small violations 
of democracy 0, 0|. For the recent status of the SM with fourth generation (SM 4 ), see e.g. 
j8 l3| and references therein. 



Although the masses of fourth generation quarks are larger than the top quark mass (the 
last analysis of the Tevatron data implies nid A > 372 GeV ll( and m U4 > 358 GeV (l2j), 



they can form bound states as a result of the smallness of the mixing between these quarks 



and ordinary SM quarks |13Hl9|. As the mass difference between these two quarks is small, 
we will refer to both members of the fourth family by m 4 . The condition for formation of 
new hadrons containing ultra- heavy quarks (Q) is given by (20| : 

. . /100 GeV \ 3/2 

For t-quark with m t = 172 GeV, Eq. (1) leads to V tq < 0.44, whereas the single top 
production at the Tevatron gives V t b > 0.74 [2lj]. When the fourth family quarks have 
sufficiently small mixing with the ordinary quarks, the hadrons made up from these quarks 
can live longer enough, and the bound state it 4 ii 4 decays through its annihilation and not 
via U4 decays to a lower family quark plus a W boson [19|. Concerning the flavor democracy 
approach, this situation is realized for parameterizations proposed in 0] and j22|, whereas 
parameterization in 0] predicts V um ~ 0.2 which does not allow formation of the fourth 
family quarkonia for m Ui > 300GeV. 

Considering the above discussions, the production of such bound states if they exist will 
be possible at LHC. The conditions for observation of the fourth SM family quarks at the 
LHC has been discussed in [l3|, 23- 3(1. As there is a possibility to observe the bound 



states which consist of fourth family quarks at the LHC, it is reasonable to investigate their 
properties, theoretically and phenomenologically. 

In the present work, we calculate the masses and decay constants of the bound state 
mesons containing two heavy quarks either both from the SM4 or one from heavy fourth 
family and the other from ordinary heavy b or c quark. Here, we consider the ground state 
mesons with different quantum numbers, namely scalars (w 4 -u 4 , u^b and u±c), pseudoscalars 
(« 4 7 5 it 4 , and S 4 7 5 c), vectors (« 4 7 M « 4 , u^j^b and u^^c) and axial vector (u^j^u^ 

-u 4 7 At 7 5 6 and UaJ^c) mesons. These mesons, similar to the ordinary hadrons, are formed 
in low energies very far from the asymptotic region. Therefore, to calculate their hadronic 
parameters such as their masses and leptonic decay constants, we need to consult some 
nonperturbative approaches. Among the nonperturbative methods, the QCD sum rules 
which is based on QCD Lagrangian and is free of model dependent parameters, is one 
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of the most applicable and predictive approaches to hadron physics. This method has been 
successfully used to calculate the masses and decay constants of mesons both in vacuum 
and at finite temperature (see for instance [32 4l|). Now, we extend the application of this 
method to calculate the masses and decay constants of the considered mesons containing 
fourth family quarkonia. The heavy quark condensates are suppressed by the inverse powers 
of the heavy quark mass. Therefore, as the first nonperturbative contributions, we take into 
account the two-gluon condensate diagrams. 

The outline of the paper is as follows. In the next section, QCD sum rules for masses 
and decay constants of the considered bound states are obtained. Section III encompasses 
our numerical analysis on the masses and decay constants of the ground state ultra heavy 
scalar, pseudoscalar, vector and axial vector mesons as well as our discussions. 



II. QCD SUM RULES FOR MASSES AND DECAY CONSTANTS OF THE 
BOUND STATES (MESONS) CONTAINING HEAVY FOURTH FAMILY 

QUARKS 

We start to this section considering sufficient correlation functions responsible for cal- 
culation of the masses and decay constants of the bound states containing heavy fourth 
generation quarks in the framework of QCD sum rules. The two point correlation function 
corresponding to the scalar (S) and pseudoscalar (PS) cases is written as: 

U s ^ = i J d 4 xe^ x (0 | T(j s(PS \x)J s ( PS \0)) | 0), (2) 

where T is the time ordering product and J s (x) = u i {x)q{x) and J PS (x) = u^x^^x) are 
the interpolating currents of the heavy scalar and pseudoscalar bound states, respectively. 
Here, the q can be either fourth family ?/ 4 quark or ordinary heavy b or c quark. Similarly, 
the correlation function for the vector (V) and axial vector (AV) is written as: 

= i J d A xe*- x (0 | T (J^ av \x) J^ AV \0)) | 0), (3) 

where, the currents = u 4: (x)'y fl q(x) and J AV = u^{x)^^q{x) are responsible for creating 
the vector and axial vector quarkonia from vacuum with the same quantum numbers as the 
interpolating currents. 

From the general philosophy of the QCD sum rules, we calculate the aforesaid correlation 
functions in two alternative ways. From the physical or phenomenological side, we calculate 
them in terms of hadronic parameters such as masses and decay constants. In QCD or 
theoretical side, they are calculated in terms of QCD degrees of freedom such as quark 
masses and gluon condensates by the help of operator product expansion (OPE) in deep 
Euclidean region. Equating these two representations of the correlation functions through 
dispersion relations, we acquire the QCD sum rules for the masses and decay constants. 
These sum rules relate the hadronic parameters to the fundamental QCD parameters. To 
suppress contribution of the higher states and continuum, Borel transformation with respect 
to the momentum squared is applied to both sides of the correlation functions. 

First, to calculate the phenomenological part, we insert a complete set of intermediate 
states having the same quantum numbers as the interpolating currents to the correlation 
functions. Performing the integral over x and isolating the ground state, we obtain 

u s ( ps) (0 I J S{PS) (0) I S(PS))(S(PS) | J s ™(0) | 0) | 

m s(PS) ~ P 2 
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where • • • represents contributions of the higher states and continuum and rris(ps) is mass of 
the heavy scalar (pseudoscalar) meson. From a similar manner, for the vector (axial vector) 
case, we obtain 

u v ( av) _ (0 | J^ AV) (0) I V(AV))(V(AV) | Jj^m I 0) ~ 
1 V - —{2 + " ' ' \°) 

rn V(AV) P 

To proceed, we need to know the matrix elements of the interpolating currents between the 
vacuum and mesonic states. These matrix elements are parametrized in terms of leptonic 
decay constants as: 



(0 | J(0) | S) = f s m s , 



(0 | J(0) | PS) = f. 



m 



2 

PS 



ps- 

m u4 + m q 

(0 | J(0) | V(AV)) = f V (AV)m v{AV) e^ (6) 

where fa are the leptonic decay constants of the considered bound state mesons. Using 
summation over polarization vectors in the V(Al^) 

" l V(AV) 

(7) 

we get, the final expressions of the physical sides of the correlation functions as: 

n s = _/>L + ... 



m| — p 2 



\2 



£2 t '"ps 
jtPS _ Ps ^m u4 ,+m q > 
2 9 
m PS ~ P 

rrV(AV) _ fv(AV) m V{AV) 

ui/ o 2 



9 If ' ^,2 



+ 



(8) 



where to calculate the mass and decay constant in the V(AV) channel, we choose the 
structure g^ v . 

In QCD side, the correlation functions are calculated in deep Euclidean region, p 2 
— Aq CD via OPE, where short or perturbative and long distance or non-perturbative effects 
are separated. For each correlation function in S(PS) case and coefficient of the selected 
structure in V(>1V) channel, we write 

n T\.p er i -\- II nonpert' (9) 

The short distance contribution (bare loop diagram in figure ([1]) part (a)) in each case is 
calculated using the perturbation theory, whereas the long distance contributions (diagrams 
shown in figure (CQ) part (b) ) are parameterized in terms of gluon condensates. To proceed, 
we write the perturbative part in terms of a dispersion integral, 

jjQCD _ J _ P_ _|_ n nonper4 , (10) 



(')) 



FIG. 1. (a): Bare loop diagram (b): Diagrams corresponding to gluon condensates. 



where, p(s) is called the spectral density. To calculate the spectral density, we calculate 
the Feynman amplitude of the bare loop diagram by the help of Cutkosky rules, where the 
quark propagators are replaced by Dirac delta function, (— 2ni)5(p 2 — m 2 ). As 

a result, the spectral density is obtained as follows: 



= _ (Jni±m,) I _ 2 m 1 + m 2 K - m 2 f 

HK ' 8vr 2V s J \ s s 2 v ' 

where + sign in (mi ± 7712) is chosen for scalar and axial vector cases and — sign is for 
psoduscalar and vector channels. Here, m\ = m U4 and m 2 is either m U4 or m c ^)- 

To obtain the non-perturbative part, we calculate the gluon condensate diagrams repre- 
sented in part (b) of figure (JTJ. For this aim, we use Fock-Schwinger gauge, x^A^x) = 0. 
In momentum space, the vacuum gluon field is expressed as: 

Al{k') = -^YG^)^\k') 1 (12) 

where k' is the gluon momentum. In the calculations, we also use the quark-gluon-quark 

vertex as: 

^ = mJ^) , (is) 
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1.1 



After straightforward but lengthy calculations, the non-perturbative part for each channel 
in momentum space is obtained as: 

™ r 1 , r 2\ Q f + e '( m i g m 2) , / U N 

^nonpert ~ J Q ) Qfa^ + ^ _ ^ _ p 2 x + ^2)4^ ^ 

where ® l (mi <H- m 2 ) means that in O l , we exchange mi and m 2 . The explicit expressions 
for 0* are given as: 

Q s = ]^x 2 \^m\x{m 2 2 {x{ll - 2x{2x{9x - 26) + 47)) + 8) 

+p 2 x{x{27x - 25) - 7){x - l) 2 ) + 2m 2 m\(m 2 2 (x{x{x(2lx - 58) + 39) 
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+ 12) - 15) - p 2 (x - l)x(x{x{7x - 13) - 3) + 12)) 
+m 2 1 (-m 2 2 p 2 (x - l)x(x(a;(2x(81x - 242) + 455) - 96) - 33) 
+mt(x(x(x(3x(36x - 145) + 652) - 414) + 72) + 15) + 3p 4 (x - l) 3 
x 2 {24x 2 - 22x - 5)) - m 2 m l {x - \){-m\p 2 [x 2 - 2)(x(14x - 27) + 15) 
+mj(3x - 5)(x(7x - 12) + 6) + p 4 (x - l)x(x(2x(7x - 13) + 3) + 12)) 
+(x - l)(-m 2 2 p 4 (x - l)x(2x(x(2x(18x - 55) + 109) - 30) - 9) 
+mjp 2 (x(x(x(x(81x - 328) + 490) - 299) + 42) + 15) 
-m\(2x - 3){x{§x{3x - 8) + 47) - 15) + 3p 6 (x - 1) V(6(x - l)x - 1)) 

+9m\(x - l) 2 x 2 {4x + 1) + 3m 2 m 5 1 x(x((8 - 7x)x + 2) - 4) I, 



-p 2 (x - l) 2 x(27x 2 - 25x - 7)) - 2m 2 m 3 (m 2 (21a: 4 - 58a; 3 + 39a; 2 + 12x - 15) 
+p 2 x(-7x 4 + 20a; 3 - 10x 2 - 15a; + 12)) + m 2 m l {x - \){m\p 2 (-14a; 4 
+27x 3 + 13x 2 - 54a: + 30) + m 4 (21x 3 - 71x 2 + 78x - 30) 
+p A x(Ux 4 - 40a; 3 + 29x 2 + 9x- 12)) + m 2 (-m 2 p 2 a;(162a; 5 - 646a; 4 + 939a: 3 
-551a; 2 + 63a: + 33) + m 4 (108a; 5 - 435a: 4 + 652a; 3 - 414a: 2 + 72a; + 15) 
+3p 4 (a; - l)V(24a; 2 - 22a; - 5)) + (a: - l)(-m 2 p 4 a;(72x 5 - 292a; 4 + 438a; 3 
-278a; 2 + 51a; + 9) + m^p 2 (81x 5 - 328a; 4 + 490a: 3 - 299a; 2 + 42a; + 15) 
+m°(-36a; 4 + 150a; 3 - 238a: 2 + 171a: - 45) + 3p 6 (a; - l) 3 x 2 (6x 2 - 6a; - 1)) 



+p 2 (x(27a; - 25) - 7)x 2 ) + 2m 2 ml(x - l) 2 x{m 2 2 (3x - 4) 

—p 2 [x — 3)x) — m 2 m\{x — l) 2 (m 2 p 2 x((7 — 2x)x — 8) + m 2 (x — l)(3x — 5) 

+p 4 x 2 (2(a; - l)x + 3)) + m\{x - l)a;(m2p 2 a;(x(-54x 2 + 56a; + 5) + 4) 

+m 4 (9(a; - l)x(4z - 1) - 8) + p 4 x 3 (24a; 2 - 22a: - 5)) + (x - l) 2 

(m 2 2 p 4 x 2 (A(7 - 6x)a; 2 + 1) + m 4 2 p 2 x(x 2 (27x - 31) - 3) 

+m 6 2 (5 - 2x(6x 2 -9x + 4)) + p 6 x 4 (6(x - l)x - 1)) 



+m 4 x(m 2 (x(17 - 2x(18(x - 3)x + 47)) + 8) + p 2 a;(x(27a; - 25) -7){x- i f) 

+m\(-m 2 2 p 2 (x - l)a;(x(a;(2a;(27x - 82) + 149) - 32) - 11) 

+m 4 2 (3x(x(x(3x(Ax - 17) + 76) - 46) + 8) + 5) + p 4 (x - l)V(24x 2 - 22a: - 5)) 

+m 2 mi(x — l)x 2 (m 2 ,p 2 (7 — a;(2x + 3)) + m 2 (3x — 5) 

+p 4 (x - l)(2(x - l)x + 3)) + (x - 1)( — m\p 4 (x — l)a;(2a;(x(2x(6x — 19) 




+9m\{x - l)V(4x + 1) + 3m 2 m^a:(7a: 3 - 8x 2 - 2x + 4) 
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+37) - 10) - 3) + m\p 2 (x{x{x{x{27x - 112) + 162) - 97) + 14) + 5) 
+ml(x(57 - 2x(3x(2x - 9) + 43)) - 15) + p 6 (x - 1) V(6(x - l)x - 1)) 

+3m 2 mlx 4 + 3mf(x - l) 2 x 2 (4x + 1) j. (15) 

The next step is to match the phenomenological and QCD sides of the correlation func- 
tions to get sum rules for the masses and decay constants of the bound states. To suppress 
contribution of the higher states and continuum, Borel transformation over p 2 as well as 
continuum subtraction are performed. As a result of this procedure, we obtain the following 
sum rules: 



J (mi+m2) 

= f S ° „ ds p PS (s) e-& + BUZ P ert, (16) 



4 rl 2 

rn^sfps 

6 M , 

(m U4 + m q ) 2 7( mi+m2 )2 



where M 2 is the Borel mass parameter and So is the continuum threshold. The sum rules 
for the masses are obtained applying derivative with respect to — to the both sides of 
the above sum rules and dividing by themselves, i.e., 



m S(PS)(V)(AV) - M a S(PSMVVAV\( \ -A , nr(S{PS){V){AV) ' ( 17 ) 



where 



i ^l+fHz^ll A* + A J (mi -h- m 2 ) . 



A s = -m 2 m 3 (a; - l)a: 2 (m 2 (14a; 2 - 29x + 14) 

+2M 2 x(7x 2 - 13a; + 6)) + m\(x - l)x 3 (m^(9x 2 - 14a: + 6) 
+3M 2 x(3x 2 - Ax + 1)) + m 2 m l (x - \){m\M 2 x 
(14a; 4 - 53a; 3 + 71a: 2 - 36a: + 6) + m 4 (7a: 4 - 28a: 3 + 40a; 2 - 25a; + 6) 
+2M 4 x 2 (Ux 4 - 40a; 3 + 29a: 2 + 9a; - 12)) + m 2 a:(m 2 M 2 x 
(-18a: 5 + 70a; 4 - 105a; 3 + 77a; 2 - 27a; + 3) + m 4 (-9x 5 + 37a; 4 
-61a: 3 + 52a; 2 - 21a; + 3) - 12M 4 a; 2 (3a: + l)(x - l) 4 ) - (x - 1) 
(-2m 2 M 4 x 3 (18a; 4 - 76a: 3 + 123a; 2 - 89a; + 24) 

+m 4 M 2 x(-9a; 5 + 40a; 4 - 71a; 3 + 68a: 2 - 33a: + 6) + m^(-3a; 5 + 14a: 4 
-27a: 3 + 29a; 2 - 15a; + 3) + 6M 6 (i - l) 3 x 3 (6x 2 - 6a; - 1)) 
-Zm\{x - l)a; 5 + m 2 mfa: 3 (7a; 2 - 8x + 1), 
A PS = -m 2 m\{x - l)x 2 (m 2 (14a: 2 - 29a: + 14) 

+2M 2 x(7x 2 - 13a: + 6)) - m\[x - l)x 3 {m 2 2 {9x 2 - 14a; + 6) 

+3M 2 x(3x 2 - 4a; + 1)) + m 2 m 1 (x - l)(m 2 M 2 x(14a: 4 - 53a; 3 + 71a: 2 - 36a: + 6) 

+m 4 (7a: 4 - 28a; 3 + 40a: 2 - 25a: + 6) + 2M 4 a; 2 (14x 4 - 40a: 3 + 29a: 2 + 9x - 12)) 
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+mjx(m 2 2 M 2 x(18x 5 - 70x 4 + 105a; 3 - 77a; 2 + 27x - 3) + mj(9x 5 - 37a: 4 
+6lx 3 - 52a; 2 + 21a; - 3) + 12MV(3a; + l)(x - l) 4 ) 
+{x - l)(-2m 2 M 4 x 3 (18a: 4 - 76a; 3 + 123a: 2 - 89a; + 24) 
+m 4 M 2 x(-9a; 5 + 40a; 4 - 71a; 3 + 68a: 2 - 33a: + 6) 
+m^(-3x 5 + 14a; 4 - 27a; 3 + 29a: 2 - 15a: + 3) 

+6M 6 (x - l)V(6a; 2 - 6x - 1)) + 3mf(x - l)x 5 + m 2 m^x 3 (7a; 2 - 8x + 1), 
A v = m 2 m\(x - l) 2 x 2 (m 2 2 (2x - 1) + 2M 2 a;(x + 2)) 

— m 4 (x — l)x 3 (m 2 (3x 2 — 3a; + 1) + M 2 (3a: — l)a: 2 ) — m 2 m 1 (a; — l) 3 a; 
(m^M 2 (2a: 2 + 3a: - 2) + m\{x - 1) + 2M 4 a:(2a: 2 - 2x + 3)) 
+m\(x - l) 2 a:(m 2 M 2 x(6x 3 - 8a; 2 + x + 2) + m^x 3 - 6x 2 + Ax - 1) 
+4M 4 x 3 (3x 2 - 2a; - 1)) + (x - l) 3 (2m 2 M 4 x 2 (-6x 3 + 10a: 2 - 3a; + 1) 
-m 4 M 2 x(3x 3 - 7a: 2 + 3x + 1) - m^(x - l) 3 + 2M 6 x 4 (6x 2 - 6a; - 1)) 
+m\x® — m 2 ml(x — l)x 4 , 
A AV = -m 2 m 3 (x - l)x 2 (m 2 (2x 2 - 5x + 2) (19) 
+2M 2 x(a; 2 - 4a; + 3)) - m\(x - l)x 3 (m 2 (3x 2 - 6x + 2) 
+M 2 x(3a: 2 - 4a; + 1)) + m 2 m 1 (x - l^tm^a^a: 3 - llx 2 + 17a; - 6) 
+m\(x 3 - Ax 2 + Ax - 1) + 2M 4 x 2 (2x 3 - 4a: 2 + 5a: - 3)) + m 2 x(m 2 M 2 x(6z 5 
-26x 4 + 43a: 3 - 31x 2 + 9x - 1) + m 4 (3a: 5 - 15a: 4 + 27x 3 - 20a; 2 + 7x-l) 
+4M 4 x 2 (3x + l)(x- l) 4 ) + (x- l)(-2m 2 M 4 z 3 (6x 4 - 28a: 3 + 45a: 2 - 31a: + 8) 
+m 4 M 2 x(-3a: 5 + 16a: 4 - 33a: 3 + 28x 2 - llx + 2) - m 6 2 (x 5 - 6x 4 + 13a: 3 - 11a: 2 
+5x - 1) + 2M 6 (s - 1) V(6x 2 - 6x - 1)) + m\(x - l)x 5 + m 2 m\{x - 1) V. 



III. NUMERICAL RESULTS 



To obtain numerical values for the masses and decay constants of the considered bound 
states containing heavy fourth family from the obtained QCD sum rules, we take the mass 
of the « 4 in the interval m U4 = (450 — 550) GeV, m b = 4.8 GeV, m c = 1.3 GeV and 
(0 | \a s G 2 | 0} = 0.012 GeV 4 . The sum rules for the masses and decay constants contain 
also two auxiliary parameters, namely Borel mass parameter M 2 and continuum threshold s§. 
The standard criteria in QCD sum rules is that the physical quantities should be independent 
of the auxiliary parameters. Therefor, we should look for working regions of these parameters 
such that our results are approximately insensitive to their variations. The working region 
for the Borel mass parameter is determined demanding that not only the higher states and 
continuum contributions are suppressed but contributions of the highest order operators 
should also be small, i.e., the sum rules for the masses and decay constants should converge. 
As a result of the above procedure, the working region for the Borel parameter is found to 
be 500 GeV 2 < M 2 < 900 GeV 2 for u A b and u 4 c, and 1200 GeV 2 < M 2 < 2000 GeV 2 
for heavy SM 4 mesons. The continuum threshold s is not completely arbitrary but 
correlated to the energy of the first exited state with the same quantum number as the 
interpolating current. We have no information about the energy of the first excitation of 
the bound states containing fourth family quarks. Hence, the only way to determine the 
working region is to choose a region such that not only the results depend weakly on this 
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parameter but the dependence of the physical observables on the Borel parameter M 2 is 
also minimal. Our numerical calculations lead to the interval (mi + m 2 + 3.3) 2 GeV 2 < s < 
(mi + m 2 + 3. 7) 2 GeV 2 for the continuum threshold. 

904 
903 

> 

O 902 



901 
900 

1200 1400 1600 1800 2000 

M 2 [GeV 2 ] 4 

FIG. 2. Dependence of mass of the scalar U4U4 on the Borel parameter, M 2 at three fixed values 
of the continuum threshold. The upper, middle and lower lines belong to the values so = (mi + 
m 2 + 3.7) 2 GeV 2 , s = (mi +m 2 + 3.5) 2 GeV 2 and s = (mi + m 2 + 3.3) 2 GeV 2 , respectively. 

904 
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> 

S 902 
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900 

1200 1400 1600 1800 2000 

M 2 [GeV 2 ] 

FIG. 3. The same as Fig. [2] but for pseudoscaler u^^Ui. 

As an example, let us consider the case of the bound state m 4 m 4 . The dependence of 
the masses of scalar U4U4, pseudoscalar W475U4, vector u^^u^ and axial vector ^4757^^4 are 
presented in figures (E115]) at three different fixed values from the considered working region 
for the continuum threshold. From these figures, we see a good stability of the masses 
with respect to the Borel mass parameter M 2 . From these figures, it is also clear that the 
results do not depend on the continuum threshold in its working region. The dependence of 
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1200 1400 1600 1800 2000 
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FIG. 4. The same as Fig. [2] but for vector 1x47^114. 
904 

903 

> 

S 902 
> 
< 

901 
900 

1200 1400 1600 1800 2000 

M 2 [GeV 2 ] 

FIG. 5. The same as Fig. [2] but for axial vector 114757^114. 

the decay constants of the scalar 114114, pseudoscalar 11475114, vector 1x47^1x4 and axial vector 
1x4757^114 are presented in figures ([MID a ^ so a ^ three different fixed values of the continuum 
threshold. These figures also depict approximately insensitivity of the results under variation 
of the Borel mass parameter in its working region. The results of decay constants also show 
very weak dependency on the continuum threshold in its working region. From a similar 
way, we analyze the mass and decay constants of the cases when one of the quarks belong 
to the heavy fourth generation and the other is ordinary bottom or charm quark. The 
numerical results deduced from the figures are collected in Tables I- VI for three different 
values of the m U4 , namely m U4 = 450 GeV, m U4 = 500 GeV and m U4 = 550 GeV. The errors 
presented in these tables are only due to the uncertainties coming from determination of 
the working regions for the auxiliary parameters. Here, we should stress that the obtained 
results in Tables I- VI are within QCD and do not include contributions coming from the 
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FIG. 6. Dependence of the decay constant of the scalar u^u^ on the Borel parameter, M 2 at three 
fixed values of the continuum threshold. The upper, middle and lower lines belong to the values 
s = (mi + m 2 + 3.7) 2 GeV 2 , s = (mi + m 2 + 3.5) 2 GeV 2 and s = (mi + m 2 + 3.3) 2 GeV 2 , 
respectively. 



5.0 



4.5 - 



O 4.0 — 



3.5 - 



3.0 



1200 



1400 



1600 
M 2 [GeV 2 



1800 



2000 



FIG. 7. The same as Fig. [6] but for the decay constant of pseudoscalar 7/4751*4. 



Higgs couplings to the ultra heavy quarks. Such contributions to the binding energy have 
been calculated in 19], where it is shown that these contributions are more than several 
GeV in the case when both quarks belong to the fourth family. The Higgs contribution 



calculated in [19] is proportional to the product of two quark masses. When we replace one 
of the ultra heavy quarks by b or c quark, the binding energy obtained in 

S3 

reduces to a 

value which is less than the QCD sum rules predictions in the present work. However, when 
both quarks belong to the fourth family, the binding energy obtained in the present work is 
very small comparing to the Higgs corrections in [l9j. 

At the end of this part, we would like to mention that the obtained QCD sum rules in the 
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FIG. 8. The same as Fig. [B]but for the decay constant of vector u^^u^. 
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FIG. 9. The same as Fig. [6] but for the decay constant of axial vector ^757^4. 



present work reproduce the masses and decay constants of the ordinary bb(cc) states when 
we set -u 4 — > b(c). The obtained numerical values in this limit are in a _good consistency with 
the existing experimental data [42| and QCD sum rules predictions (iol . 41 . 



To sum up, against the top quark, the heavy fourth generation of quarks that have 
sufficiently small mixing with the three known SM families form hadrons. Considering the 
arguments mentioned in the text, the production of such bound states will be possible at 
LHC. Hoping for this possibility, we calculated the masses and decay constants of the bound 
state objects containing two quarks either both from the SM 4 or one from heavy fourth 
generation and the other from observed SM bottom or charm quarks in the framework of 
the QCD sum rules. The obtained numerical results approach to the known masses and 
decay constants of the bb and cc heavy quarkonia, when the fourth family quark is replaced 
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mass (GeV) 


U4C 


U46 


U4U4 


Scalar 


453.01 ±0.25 


456.45 ± 0.25 


901.68 ±0.50 


Pseudoscalar 


452.62 ±0.15 


455.95 ±0.15 


901.12 ±0.30 


axial vector 


453.00 ± 0.25 


456.44 ± 0.25 


901.70 ±0.50 


vector 


452.62 ±0.15 


455.94 ±0.15 


901.13 ±0.30 



TABLE I. The values of masses of different bound states obtained using m Ui = 450 GeV. 



mass (GeV) 


U4C 






Scalar 


502.91 ± 0.28 


506.36 ±0.28 


1001.61 ±0.55 


Pseudoscalar 


502.52 ±0.17 


505.86 ±0.17 


1001.04 ±0.33 


Axial Vector 


502.91 ± 0.28 


506.35 ±0.28 


1001.60 ±0.55 


Vector 


502.57 ±0.17 


505.85 ±0.17 


1001.04 ±0.33 



TABLE II. The values of masses of different bound states obtained using m Ui = 500 GeV . 
by the bottom or charm quark. 
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Leptonic decay constant f (GeV) 


U4C 


u^b 


U4U4 


Scalar 


0.10 ±0.01 


0.12 ±0.01 


0.26 ± 0.03 


Pseudoscalar 


0.14 ±0.01 


0.27 ±0.01 


4.19 ±0.20 


Axial Vector 


0.10 ±0.01 


0.12 ±0.01 


0.26 ± 0.03 


Vector 


0.14 ±0.01 


0.27 ±0.01 


4.18 ±0.20 



TABLE VI. The values of decay constants of different bound states obtained using m U4 = 550 GeV. 
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